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ABSTRACT 

This paper explores how orbits in a galactic potential can be impacted by large am¬ 
plitude time-dependences of the form that one might associate with galaxy or halo 
formation or strong encounters between pairs of galaxies. A period of time-dependence 
with a strong, possibly damped, oscillatory component can give rise to large amounts 
of transient chaos, and it is argued that chaotic phase mixing associated with this tran¬ 
sient chaos could play a major role in accounting for the speed and efficiency of violent 
relaxation. Analysis of simple toy models involving time-dependent perturbations of 
an integrable Plummer potential indicates that this chaos results from a broad, pos¬ 
sibly generic, resonance between the frequencies of the orbits and harmonics thereof 
and the frequencies of the time-dependent perturbation. Numerical computations of 
orbits in potentials exhibiting damped oscillations suggest that, within a period of 
10 dynamical times to or so, one could achieve simultaneously both ‘near-complete’ 
chaotic phase mixing and a nearly time-independent, integrable end state. 
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1 MOTIVATION 

Dating back at least to the pioneering work of Lynden-Bell 
(1967), it has been recognised that violent relaxation, viewed 
broadly as the mechanism whereby, even in the absence 
of dissipation, a galaxy evolves towards an equilibrium or 
near-equilibrium state, must be a collective process some¬ 
how involving phase mixing. However, not just any phase 
mixing will do: the simplest time-independent models which 
one might envision, such as Lynden-Bell’s balls rolling in 
a pig-trough, do not approach an equilibrium nearly fast 
enough (c/. Kandrup 2001). For this reason, it has often 
been assumed that the fact that the potential is strongly 
time-dependent is crucial in making the mixing sufficiently 
efficient. 

In the mid-1990’s, however, it was recognised that it 
may not be time dependence per se that is crucial. Rather, 
what seems essential is that the flow be chaotic, exhibit¬ 
ing an exponentially sensitive dependence on initial condi¬ 
tions. The important point here (c/. Kandrup and Mahon 
1994, Mahon, Abernathy, Bradley, and Kandrup 1995, Mer- 
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ritt and Valluri 1996) is that, whereas the phase mixing of 
regular orbits proceeds as a modest power law in time, the 
phase mixing of chaotic orbits proceeds exponentially. Thus, 
e.g., a localised ensemble of initial conditions corresponding 
to chaotic orbits will, when evolved in a time-independent 
Hamiltonian, begin by diverging exponentially from its lo¬ 
calised state and then converge exponentially towards an 
equilibrium or near-equilibrium state (Kandrup 1998). 

Unfortunately, however, most work on chaotic phase 
mixing hitherto has focused on flows in time-independent 
Hamiltonian systems, possibly allowing for low amplitude 
perturbations (c/. Kandrup, Pogorelov, and Sideris 2000) 
but ignoring the possibility of a large amplitude time- 
dependence. This work on (nearly) time-independent sys¬ 
tems confirms that chaotic phase mixing can be very effi¬ 
cient. However, it is clear that this phenomenon can be im¬ 
portant only if the potential admits large measures of chaotic 
orbits. The obvious point, then, is that time-independent 
density distributions corresponding to galactic equilibria 
or near-equilibria are not likely to contain huge measures 
of chaos. As stressed, e.g., by Binney (1978), interestingly 
shaped equilibria which might be expected to admit global 
stochasticity must contain large measures of regular orbits 
to serve as the ‘skeleton’ of the interesting structure. 

The obvious questions, then are: might one expect 
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to have much larger measures of chaos in realistic time- 
dependent systems? and should these orbits exhibit the qual¬ 
itative behaviour which one would expect to observe during 
an epoch of violent relaxation? The aim of this paper is to 
argue that the answers to these questions is: Yes! 

The crucially important point is that allowing for a 
time-dependent Hamiltonian can lead to resonances that 
trigger an epoch of transient chaos: even if the initial and fi¬ 
nal density distributions are nearly integrable and admit no 
global stochasticity, the intervening states could still corre¬ 
spond to orbits exhibiting significant exponential sensitivity. 

As a paradigmatic example, one can consider a test par¬ 
ticle moving in a spherically symmetric, constant density 
mass distribution which is, in addition, subjected to a pe¬ 
riodic sinusoidal perturbation corresponding to a force per 
unit mass F = —f3r cosuit. In the absence of the perturba¬ 
tion, a particle will feel an isotropic harmonic oscillator po¬ 
tential and, as such, the separation 5r between two nearby 
orbits will execute simple harmonic motion. However, the 
presence of the periodic driving complicates matters, lead¬ 
ing to a more complex equation for Sr, a typical component 
of which will satisfy 

(ESr o o 

, - + Q,^5r = 0, with Or = a + jScosut, (1) 

dt^ 

a classic example of a Hill equation. The obvious point, then, 
is that even if a is positive and |/3| < a, Sr can grow exponen¬ 
tially in time, corresponding to a chaotic orbit. The periodic 
driving has triggered a parametric instability. 

Real systems are not constant density, so that the un¬ 
perturbed (time-independent) orbits are more complex, and 
the time-dependent perturbations to which the system is 
subjected may not be strictly sinusoidal. However, as will 
be shown in this paper, one might still expect generically 
that a system exhibiting damped oscillations of the form as¬ 
sociated with many scenaria involving collective relaxation 
will exhibit a large amount of chaos. 

That time-dependent perturbations can lead to tran¬ 
sient chaos has been recognised already in other branches of 
physics, including nonneutral plasmas (cf. Qian, Davidson, 
and Chen 1999, Strasburg and Davidson 2000) and charged 
particle beams (Gluckstern 1994). Indeed, Gluckstern has 
suggested that chaos triggered by time-dependent perturba¬ 
tions may be responsible for undesirable broadening (‘emit- 
tance growth’) in a focused charged particle beam by eject¬ 
ing particles from the central core into an outerlying halo. 

Section 2 of this paper begins by discussing some ‘natu¬ 
ral’ physical expectations one might have regarding the pos¬ 
sibility of transient chaos induced by time-dependent pertur¬ 
bations, and then describes a set of numerical experiments 
that were performed to test these expectations. Section 3 fo¬ 
cuses on the possibility of transient chaos in its purest form, 
considering in detail what can happen when an integrable 
spherical system is subjected to an undamped oscillatory 
perturbation with a single frequency. Section 4 extrapolates 
from this example to consider more realistic computations 
that capture two important aspects associated with violent 
relaxation: (1) Allowing for damped oseillations. Can one 
trigger efficient chaotic mixing within (say) a period of 10 
dynamical times to or so while simultaneously damping the 
original strongly time-dependent potential to a nearly time- 
independent state? (2) Allowing for a variable pulsation fre¬ 


quency. As the density distribution changes, one might ex¬ 
pect that the pulsation frequency will also change, and the 
obvious question is: do such changes matter significantly? 
Section 5 summarises the principal conclusions and describes 
extensions of the work described herein currently underway. 

2 WHAT WAS EXPECTED AND WHAT WAS 

DONE 

2.1 Physical expectations 

As stated already, the basic idea is that an appropriate time- 
dependence added to a potential can trigger chaos via a 
parametric resonance: if orbits are subjected to a perturba¬ 
tion with power at appropriate frequencies, one might expect 
otherwise regular orbits to become chaotic. In particular, if 
unperturbed orbits have power concentrated at frequencies 
~ II and the perturbation has significant power at frequen¬ 
cies u sufficiently close to Q, there is the possibility of reso¬ 
nance overlap which, as is well known {cf. Lichtenberg and 
Lieberman 1992), can lead to the onset of global stochastic¬ 
ity. 

As a particularly simple example, one knows that if a 
regular orbit with substantial power at some frequency 12 
is perturbed by a sinusoidal perturbation with frequency 
u = Q, the orbit will typically exhibit a drastic response 
which can manifest sensitive dependence on initial condi¬ 
tions. This is, e.g., the physics responsible for resonant cir¬ 
cuits in elementary electronics. The real issue, one might 
argue, is simply: ‘how close’ must 12 and u be to trigger 
strong exponential sensitivity? 

Detailed analysis of the effects of resonance overlap sug¬ 
gests generically that the width of such resonances will be an 
increasing function of amplitude. Even if a slightly ‘detuned’ 
small amplitude perturbation has no apparent effect, a large 
perturbation with identical frequency may elicit a huge re¬ 
sponse. Indeed, the numerical models described in Sections 

3 and 4 of this paper indicate that, for large amplitude per¬ 
turbations - fractional amplitude of order 10 — 20% or more 
- the resonance can be very broad. In particular, these mod¬ 
els reveal that, if the power in the orbits is concentrated at 
frequencies ~ 12, one can get a nontrivial increase in chaos - 
in both the relative number of chaotic orbits and the value 
of the largest Lyapunov exponent - for sinusoidal perturba¬ 
tions with 0.112 u 3012. In other words, the resonance 
can be more than two orders of magnitude wide! 

The important point, then, is that if resonances of this 
sort are really so broad, one might expect transient chaos to 
be extremely common, if not ubiquitous, in time-dependent 
galactic potentials. Numerical simulations of galaxy encoun¬ 
ters and mergers and most models of galaxy and halo for¬ 
mation imply that a system approaching equilibrium will 
exhibit damped oscillations. To the extent, however, that 
one is dealing with collisionless relaxation, there is dimen¬ 
sionally only one natural time scale in the problem, namely 
the dynamical i_D ~ If^/Gp, with p a characteristic density. 
In particular, to should set the oscillation time scale as well 
as the orbital time scale. The exact numerical values of these 
time scales will involve numerical coefficients which will in 
general be unequal and vary as a function of location within 
the galaxy. If, however, one only needs to assume that the 
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oscillation and orbital time scale agree to within an order of 
magnitude or so, it would seem likely that the conditions ap¬ 
propriate for resonance-induced transient chaos could arise 
in much, if not all, of the galaxy. 

The crucial recognition, then, is that even a relatively 
short period of exponential sensitivity can result in com¬ 
paratively efficient chaotic phase mixing. Earlier analysis of 
flows in time-independent potentials (c/. Mahon, Abernathy, 
Bradley, & Kandrup 1995, Merritt & Valluri 1996, Kandrup 
1998) have demonstrated that the presence of chaos can dra¬ 
matically enhance the degree to which a system becomes 
‘shuffled’ and, hence, the rate at which a localised orbit en¬ 
semble evolves towards an equilbrium. By complete analogy, 
one might expect that transient chaos in a self-consistently 
determined time-dependent potential will again enhance the 
efficacy of ‘mixing’ and, by so doing, facilitate a more rapid 
approach towards an equilibrium. 

It should be noted explicitly that the physical mecha¬ 
nism for the generation of chaos described here is very dif¬ 
ferent from the mechanism whereby the amount of chaos 
can increase - or decrease - in more slowly varying poten¬ 
tials (c/. Kandrup and Drury 1998, Contopoulos 2002). If 
the potential is slowly varying, one can visualise an orbit as 
moving in a phase space which is slowly changing but which 
is nearly constant over a time scale ~ to. It is then natu¬ 
ral to suppose that at some times the orbit finds itself in a 
phase space region that is ‘chaotic’ in the sense that there 
is a sensitive dependence on initial conditions, but that at 
other times it finds itself in ‘regular’ regions where there is 
no such sensitive dependence. In the setting envisioned in 
this paper, the phase space is changing on a time scale com¬ 
parable to - or even shorter than - the orbital time scale, so 
that this picture is necessarily lost. Moreover, in the setting 
described here one would expect generically a systematic 
increase in the amount and degree of chaos, not simply a 
change which could be either positive or negative. 

Finally, it should be noted that, to a certain extent, the 
term transient chaos is necessarily ‘fuzzy.’ Chaos, like ordi¬ 
nary Lyapunov exponents, is only defined in an asymptotic 
t ^ oo limit. The point, however, is that, even over finite 
time intervals, it is physically meaningful to ask whether or¬ 
bits exhibit an exponentially sensitive dependence on initial 
conditions. This is, e.g., the notion that motivates the def¬ 
inition of finite time Lyapunov exponents (c/. Grassberger, 
Badii, and Poloti 1988), which have become accepted tools in 
nonlinear dynamics. From a phenomenological point of view, 
it makes perfect sense to identify an orbit as exhibiting sig¬ 
nificant transient chaos if, for some finite time > t d , an orbit 
exhibits significant exponential sensitivity. In many cases, it 
is clear by visual inspection and/or through a computation 
of a finite time Lyapunov exponent whether such transient 
chaos is present. In other cases, however, things are not so 
completely clear. For example, for the models considered in 
this paper it proves difficult to determine with complete reli¬ 
ability a minimum frequency at which transient chaos ‘turns 
on’ since this ‘turning on’ is comparatively gradual. 


2.2 Numerical experiments 

The numerical computations described in this paper in¬ 
volved perturbations of a Plummer sphere, which is char¬ 
acterised (in units with G = m — 1) hy the potential 


Vo{x,y,z) = - 


{1 + x '^ -\- + 2 ; 2 ) i /2 


( 2 ) 


This form was selected (1) for its inate simplicity and, es¬ 
pecially, (2) because orbits in this potential are all strictly 
integrable. Any chaos that is observed can be associated un¬ 
ambiguously with the perturbations that were introduced. 
These were assumed to take the form 


Vi{x,y,z,t) = - 


(1 -I- -I- -I- ’ 


( 3 ) 


where the quantities m and a are allowed to vary in time. 
A variety of different time-dependences were considered, in¬ 
cluding: 

• strictly sinusoidal oscillations in m or o, i.e., 


m = mo sin tot and a = ao + 5a sin cot. 


( 4 ) 


This allowed one to focus on the physical mechanism in its 
‘purest’ form. 

• sinusoidal oscillations that damp exponentially or as a 
power law in time, i.e., 

m — moe~“* sin cut (5) 


or 


m — mo 


sin cut 
{to + t)P ’ 


( 6 ) 


with p = 1 or 2, and analogous formulae for a{t). This al¬ 
lowed one to confirm that the physical effects associated with 
a purely oscillatory perturbation persist if, as is natural in 
the context of violent relaxation, one allows for a galaxy that 
damps towards equilibrium. 

• nonoscillatory damping, i.e., 


— at 

m = moe 


( 7 ) 


mo 


{to + t)P 


( 8 ) 


and analogous formulae for a{t). This allowed one to demon¬ 
strate that the presence of an oscillatory component is cru¬ 
cial in order to get a substantial amount of transient chaos. 
As will be described below, computations with perturba¬ 
tions given by eq. (7) or (8) typically exhibited compara¬ 
tively minimal amounts of transient chaos. 

• variable driving frequency, assuming that 


cu(t) = cuo-I-5(t), (9) 

with 5{t) a randomly varying variable that samples an 
Ornstein-Uhlenbeck process {cf. Van Kampen 1981). What 
this means is that 5 is treated as Gaussian coloured noise, 
so that its statistical properties are determined uniquely by 
its first two moments, which are assumed to satisfy 

{5{t)}=0 and {5{ti)5{t2)) ^ A'^ exp{—\ti — t 2 \/tc). (10) 


Here A represents the typical ‘size’ of the random compo¬ 
nent and tc the time scale over which it changes apprecia¬ 
bly. Considering perturbations of this form allowed one to 
confirm that, because the resonance is comparatively broad, 
permitting the driving frequency to drift (within limits) has 
a relatively minor - albeit potentially significant - effect. 
This again is important physically. Since the density distri¬ 
bution varies as a galaxy evolves towards an equilibrium. 
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one would anticipate that the characteristic oscillation fre¬ 
quency (or frequencies) will change. 

Computations were performed for ensembles of 1600 ini¬ 
tial conditions which were generated in two different ways. 
(1) Broad ‘representative’ ensembles were generated by uni¬ 
formly sampling a constant energy hypersurface. A study 
of orbits generated from such ensembles allowed one to de¬ 
rive systematic effects which might be expected to act in a 
galaxy as a whole. (2) Localised ensembles were generated 
by sampling small phase space hypercubes. Studying orbits 
generated from such ensembles allowed one to determine the 
extent to which the ‘smooth’ behaviour associated with the 
representative ensembles reflected the choice of ensemble as 
opposed to properties of the resonance. 

Most experiments involved perturbations at the 50% 
level or less, i.e., mo 0.5 and/or Sa ^ 0.5. For the in¬ 
termediate energies considered throughout, a ‘typical’ or¬ 
bital frequency (derived from a Fourier transform) was 
n ~ 0.3 — 0.35, corresponding to to = 2ti/Q. ~ 20 in ab¬ 
solute units. On physical grounds one might expect that 
the perturbations should damp on a time scale compara¬ 
ble to, but somewhat longer than, to, so that most exper¬ 
iments assumed that (to within factors of a few) to and 
a~^ ~ 5tD ~ 100. (Computations of chaotic phase mixing in 
time-independent potentials [cf. Kandrup & Novotny 2002] 
typically exhibit an exponential approach towards a near¬ 
invariant distribution on a time scale ~ 5 to!) 


3 PARAMETRIC RESONANCE AND 
TRANSIENT CHAOS 

Computations involving representative orbit ensembles have 
led to several unambiguous conclusions: 

1. Integrations which do not involve large scale time vari¬ 
ations on a time scale ~ to yield little, if any, transient 
chaos. In particular, as will be illustrated more carefully in 
the following Section, allowing for nonoscillatory damped 
perturbations of the form given by eq. (7) do not result in 
large measures of chaos, even for large amplitudes mo ^ 0.5. 
However, allowing for oscillations in m and/or a can trigger 
substantial amounts of chaos. 

2. There are at least two solid reasons to believe that this 
chaos is triggered by a resonance. More obvious, perhaps, is 
the fact that the largest increases in chaos - as reflected by 
both the number of chaotic orbits and the size of a typical 
Lyapunov exponent - arise for driving frequencies compara¬ 
ble to, or somewhat larger than, the frequencies for which 
the unperturbed orbits have most of their power. Equally 
significant, however, is the fact that the driving frequencies 
which result in the greatest amount of chaos correspond pre¬ 
cisely to those frequencies for which the orbital energies are 
‘shuffled’ the most, i.e., where different particles with the 
same initial energy end up with the largest spread in final 
energies. A priori there need be no direct connection be¬ 
tween increases in chaos and large-scale shuffling of energies. 
However, one would expect resonant couplings to lead to sig¬ 
nificant changes in energies. A correlation between changes 
in energy and the amount of chaos thus corroborates the 
interpretation that this chaos is resonant in origin. 

3. At least for the case of large amplitude oscillations, 
mo ^0.1 or so, the resonance is comparatively broad. One 




Figure 1. (a) The relative fraction / of chaotic orbits in a repre¬ 
sentative 1600 orbit ensemble which are subjected to strictly sinu¬ 
soidal oscillations with variable frequency ui. The different curves 
represent, from bottom to top, driving amplitudes mo =0.1, 0.2, 
0.3, 0.4, and 0.5. (b) The same data plotted as a function of 
logio iO. 


observes significant increases in both the relative measure 
of chaotic orbits and the size of the largest Lyapunov ex¬ 
ponent whenever the driving frequency lo and the natural 
orbital frequencies Q agree to within an order of magnitude 
or so. If, e.g., one considers an energy where the natural fre¬ 
quencies peak for Qrn ~ 0.3 — 0.35 or so, noticeable increases 
in the degree of chaos are observed for 0.035 ^ to ^ 10.0, i.e., 
0.1 ~ Uj/D.m ~ 30. 

4. The width of the resonance, i.e., the range of values of to 
for which one sees an appreciable increase in the degree of 
chaos, appears to be an increasing function of amplitude. 

5. The resonance appears to be smooth in the sense that, 
for fixed amplitude, the relative number of chaotic orbits 
and the size of the largest Lyapunov exponent both vary 
smoothly as functions of the driving frequency u>. In partic¬ 
ular, plots of the fraction / of the orbits that are chaotic or 
the mean size (x) of the largest Lyapunov exponent as func¬ 
tions of driving frequency typically exhibit a single peak. 
There is no obvious evidence for complex structure within 
the resonance. 

6. For fixed to, the relative number of chaotic orbits and the 
size of the largest Lyapunov exponent are both an increas¬ 
ing function of perturbation amplitude, although it appears 
that the fraction / of orbits that are chaotic can asymptote 
towards a fixed value /m for sufficiently large amplitude. 
For some values of w, especially those near the center of the 
resonance, this /m can approach unity; for other choices of 
frequency, /m can be substantially smaller. 

7. Whether or not the perturbation breaks spherical sym¬ 
metry appears to be relatively unimportant, e.g., making 
the parameter a time-dependent or otherwise different from 
unity does not change things all that much. Provided that 
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Figure 2. (a) The mean value (x) of the chaotic orbits in the 
ensembles used to generate the preceding Figure, again plotted 
as a function of to. (b) The same data plotted as a function of 
logio UJ. 

the system manifests large amplitude oscillations, one can 
apparently get just as much chaos for spherical systems as 
for nonspherical systems. 

Evidence for several of these conclusions is provided in 
Figures 1-4, all generated from the same set of 1600 initial 
conditions, evolved for a total time t = 512 in the presence of 
a strictly sinusoidal perturbation of the form given by eq. (4) 
with 0 = 1. The two panels of Figure 1 exhibit the relative 
measure / of chaotic orbits as a function of driving frequency 
u for five different choices of amplitude mo. The top panel 
exhibits /(to) on a linear scale, thus allowing one to focus 
on the behaviour of the resonance at higher frequencies; the 
lower exhibits / as a function of log ,q to, which allows one 
to see more clearly the behaviour at lower frequencies. 

Most obvious, perhaps, from Figure 1 is the fact that 
the width of the resonance is an increasing function of mo, 
although this increase appears to be comparatively small 
for mo > 0.2 or so. Also evident is the fact that, except 
perhaps for the lowest amplitude, mo = 0.1, the fraction 
/ appears to vary smoothly with frequency. In each case, / 
peaks at a frequency ujmax comparable to the typical natural 
frequencies ~ 0.3 — 0.35 associated with the unperturbed 
orbits. However, oJmax does seem to vary somewhat as a func¬ 
tion of mo. For larger amplitudes, near uimax essentially all 
the orbits are chaotic, i.e., fiumax) ~ 1. For lower frequen¬ 
cies fioJmax) can be significantly less than unity. Overall, 
there is evidence for significant amounts of transient chaos 
for 0.35 ^ cu ^ 10.0. 

The two panels of Figure 2 exhibit (x), the mean value 
of the largest finite time Lyapunov exponent for the same en¬ 
sembles, again plotted on both linear and logarithmic scales. 
Here (x) was extracted by first identifying those orbits in the 
ensemble that were deemed chaotic and then computing the 
mean value of x for those orbits. 


Figure 3 exhibits SE, the root mean spread in energies 
at time t = 512 for the chaotic orbits identified in Figures 
1 and 2. It is clear visually that, as one would expect if the 
transient chaos is associated with a resonant coupling, the 
frequencies which result in the largest measures of chaotic 
orbits and the largest (x) also result in the largest shifts in 
energy. However, it is also evident that, especially for the 
higher amplitude perturbations, SE exhibits a more com¬ 
plex dependence on u than do /(w) or (x(w)). This does 
not contradict the assertion that transient chaos is induced 
by a resonant coupling, but it does suggest that the amount 
and degree of chaos in the broad resonance region is less 
sensitive to the pulsation frequency than is the shuffling in 
energies. This has potentially significant implications for vi¬ 
olent relaxation, where one is interested in both (i) shuffling 
the energies of the individual masses and (ii) randomising 
their phase space locations on a constant energy hypersur¬ 
face. 

One other point is also evident from this Figure, namely 
that, at least for frequencies well into the resonance, the 
spread in energies can be very large. This reflects the fact 
that, because of the resonant coupling, a significant fraction 
of the orbits have acquired large energies that set them along 
trajectories involving excursions to very large radii, r > 10 
or more. Allowing for a resonance that is this strong to last 
for as long as t ~ 20tD is completely unrealistic in the con¬ 
text of any model for violent relaxation. The computations 
described in Section 4b, which incorporate strong damping, 
avoid this artificial behaviour. 

In part because of this effect, any uniform prescrip¬ 
tion used to distinguish between ‘regular’ and ‘chaotic’ or¬ 
bits is necessarily somewhat ad hoc and, as such, fraught 
with difficulties. Figs. 1 and 2 reflect an analysis in which 
‘chaotic’ was defined as corresponding to a finite time Lya¬ 
punov exponent for the interval 0 < f < 512 with a value 
X > 0.012, arguably a fairly large threshhold value. Lower¬ 
ing the threshhold value of 0.012 to a smaller value would 
of course increase the estimated fraction / of chaotic orbits 
and decrease the mean value of (x) for those orbits deemed 
chaotic. In any event, the choice x > 0.012 implies that there 
is little if any chaos for ui ^ 0.01 and w ^ 10, which is con¬ 
sistent with the fact that localised ensembles exhibit little 
if any evidence of chaotic phase mixing for w < 0.01 and 
w > 10 or so. Despite this, however, it should be recognised 
explicitly that this prescription really identifies orbits which 
are ‘obviously’ or ‘strongly’ chaotic, and that it may be ig¬ 
noring a considerable number of orbits which are only very 
weakly chaotic. 

This prescription, like any other which involves a fixed 
minimum value and a fixed integration time, is open to crit¬ 
icism. For those frequencies and amplitudes where the res¬ 
onance is especially strong and many orbits quickly achieve 
large energies, it might seem more appropriate to consider a 
shorter time interval, during which the orbits are restricted 
to smaller radii. However, for frequencies and amplitudes 
for which the resonance is weaker, one might wish instead 
to integrate for much longer times so as (hopefully) to allow 
for clear distinctions to be made between regular and ‘ex¬ 
tremely sticky’, nearly regular (c/. Contopoulos 1971) orbits. 
The choice of a time interval t = 512 represents a compro¬ 
mise necessitated by the idealised nature of the computa¬ 
tions described in this section. 
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Figure 3. (a) The root mean squared spread in energies, 
SErms(t = 512), for the orbit ensembles used to generate the first 
two Figures, (b) The same data plotted as a function of logj^g oj. 
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Figure 4. The relative fraction / of chaotic orbits in a representa¬ 
tive 1600 orbit ensemble which are subjected to strictly sinusoidal 
oscillations of the form given by eq. (4) of variable amplitude mg 
with u) = 1.4 (upper curve), u> = 3.5 (middle curve), and u> = 7.0 
(lower curve). 
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In any event, what is inevitable is that, for chaotic or¬ 
bits in any potential, one must expect correlations between 
the energy of the orbit and the size of a typical finite time 
Lyapunov exponent x- For the potential considered here, 
one finds typically that the orbits deemed regular are pre¬ 
cisely those which exhibit the smallest changes in energy 
and hence, overall, the smallest excursions from the origin. 
However, for those orbits which are chaotic, those spending 
more time at larger radii and which, for this reason, have 
larger orbital time scales, tend to have somewhat smaller 
values of x- 


Figure 4 exhibits the relative measure of chaotic orbits 
as a function of amplitude mo for three different frequen¬ 
cies, namely w = 1.4, 3.5, and 7.0. For w = 1.4, a frequency 
in the middle of the resonance, for amplitudes as large as 
mo = 0.25 almost all the orbits exhibit evidence of chaos. 
For u! = 3.5, a value somewhat closer to the edge of the 
resonance, the relative abundance of chaotic orbits is some¬ 
what smaller. However, for mo = 0.25 as many as 75% of 
the orbits are clearly chaotic, a fraction which does not in¬ 
crease significantly if mo is increased. For to = 7.0, a value 
near the edge of the resonance, the fraction of the orbits 
which is chaotic is smaller and does not appear to level off 
for mo > 0.25 or so. Indeed, if the amplitude is increased to 
a value as large as mo = 1.0 as many as 50% of the orbits 
prove chaotic, i.e., /(mo = 1) « 0.5. 

The uniform criterion used to identify chaotic orbits for 
Figures 1 and 2 fails for mo 0.1. For that reason. Figure 
4 implemented a different, somewhat more subjective crite¬ 
rion which has been used successfully in distinguishing be¬ 
tween regular and ‘very sticky’ chaotic orbits in the triaxial 
generalisations of the Dehnen potentials (Siopis & Kandrup 
2000). What this entailed was ordering the 1600 computed 
values of x a-t various times Ta, plotting the ordered list 
of x{Ta)’s for different values of Ta, and then in each case 
searching for a ‘kink’ in the curve. This prescription appears 
to permit an accurate determination of a fraction / of the 
orbits exhibiting chaotic behaviour which is insensitive to 
the precise choice of Ta, although the estimated size of the 
finite time x{Ta) does depend on the sampling interval. 

Additional insights into the nature of the resonance can 
be derived from Figure 5, which shows the ^-components of 
the composite (normalised) Fourier transforms of coordinate 
and force per unit mass, |a:(II)| and |aa;(fl)|, constructed 
by combining spectra for 1600 unperturbed w = 0 orbits 
generated from the same initial conditions used to gener¬ 
ate Figures 1-3. (Because of spherical symmetry the y- 
and 2 -components are statistically identical.) As noted al¬ 
ready, for this particular choice of energy, |a;(fl)| peaks at 

~ 0.3 — 0.35 and is comparatively negligible for much 
higher and lower frequencies. By contrast, |aa;(fl)| has multi¬ 
ple peaks associated with higher frequency harmonics which 
arise because the force per unit mass is a nonlinear function 
of position. It is natural to suppose that it is these harmonics 
that are responsible for the comparatively large responses at 
frequencies II large compared with Q.m that are evident in 
Figures 1 - 3. It is particularly interesting that this response 
cuts off for driving frequencies ui large compared with 20, 
rather than O, a fact that would suggest that one is observ¬ 
ing the effects of a broad 2 : 1 resonance similar to that 
which arises in the standard Matthieu equation (see, e.g., 
Matthews and Walker 1964). 

Much of the smoothness in plots such as those given 
in Figure 2 reflects the fact that the data were generated 
from ‘representative’ ensembles of initial conditions. If the 
same analysis is repeated for different localised ensembles of 
initial conditions, even ensembles with the same energy, one 
observes signficant variability in both the size of a typical fi¬ 
nite time Lyapunov exponent and the value of the frequency 
uj for which (x) is maximised. And similarly, plots of quan¬ 
tities like (x) as a function of ui for such ensembles exhibit 
more structure than do plots for representative ensembles. 
This presumably reflects the fact that the natural frequen¬ 
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Figure 5. (a) The power spectrum |a:(r2)| associated with the x- 
component of the orbits generated from the 1600 initial conditions 
used to created Figs. 1-3, evolved in the absence of any time- 
dependent perturbation and plotted as a function of 2Q. (b) The 
power spectrum |aa;(f2)| associated with the x-component of the 
force per unit mass associated with the same orbits. 

cies of the otherwise regular orbits vary with phase space 
location, and that some frequencies are more susceptible to 
the resonance than others. This additional structure could 
have important implications for near-equilibria subjected to 
nearly periodic perturbations, which could impact certain 
orbits much more than others. However, one might expect 
that, when considering the ‘global’ properties of violent re¬ 
laxation, such details would tend to wash out. 


4 A TOY MODEL FOR VIOLENT 
RELAXATION 

4.1 Variable pulsation frequencies 

Allowing for a randomly varying frequency of the form given 
by eqs. (9) and (10) leads to several significant conclusions. 

1. Perhaps the most important is that allowing for a compar¬ 
atively slow frequency drift can actually increase the relative 
measure of chaotic orbits. Suppose that the unperturbed fre¬ 
quency cJo is well into the resonant region and that A is not 
sufficiently large to push u> outside. In this case, the intro¬ 
duction of a perturbation h(t) with an autocorrelation time 
tc long compared with tn tends to convert many, if not all, 
the ‘regular’ orbits into orbits with appreciable exponential 
sensitivity. This behaviour can be understood by supposing 
that, in the absence of the frequency drift, there exist some 
regular orbits which, because of the particular values of their 
natural frequencies, barely manage to avoid resonating with 
the perturbation. Allowing for a drift facilitates an improved 
resonance coupling which can make (some of) these orbits 
chaotic for at least some of the time. 


Figure 6. (a) Finite time Lyapunov exponents for 1600 represen¬ 
tative orbits evolved for a time t = 512 in an undamped Plum¬ 
mer potential pulsed with amplitude m = 0.1 and a frequency 
frequency uj(t) = uiq + ^(t), with cjq = 3.5, autocorrelation time 
tc = 320 = 16t£), and = 0.35 (dashed curve) and = 2.8 
(triple-dot-dashed). The solid curve contrasts exponents derived 
for an evolution with (5 = 0. (b) The same for tc = 160. (c) tc = 80. 
(d) tc = 40. (e) tc = 20. (f) tc = 10. 


2. Provided that cJo and wo ± A are well within the reso¬ 
nance region, it is also apparant that, as the autocorrelation 
time tc decreases, the degree of chaos, as probed by the mean 
(x), tends to decrease, even though the relative fraction of 
chaotic orbits may well remain equal to unity. 

3. For fixed autocorrelation time tc, the mean (x) mani¬ 
fests only a relatively weak dependence on A, at least for 
0.1 A/cJo 0.8. For some choices of tc increasing A causes 
an increase in the mean finite time x; in other cases increas¬ 
ing A decreases the mean. However, the observed variations 
are invariably small. 

All three of these points are illustrated in Figures 6 and 
7. The dashed and triple-dot-dashed curves in Figure 6 cor¬ 
respond to plots of ordered hnite time Lyapunov exponents 
X associated with 1600 initial conditions subjected to un¬ 
damped pulsations with variable frequency cj = cuo -b 5(t) 
for mo = 0.1, cjJo = 3.5 and different choices of A and tc- 
The solid curve in each panel contrasts the exponents gen¬ 
erated for the same initial conditions evolved with (5 = 0. It 
is evident that, for the largest values of tc, the differences 
between the curves with h zero and nonzero are relatively 
small, at least for the largest values of X- However, it is also 
apparent that the time-dependent frequency drift has sig- 
nihcantly reduced the relative number of orbits with very 
small finite time x’s- Also evident is the fact that decreasing 
the autocorrelation time tc, i.e., making the frequency drift 
more quickly, tends generically to decrease the typical value 
of Xi although the relative measure of chaotic orbits may 
not change appreciably. For this particular choice of cjJq and 
mo, increasing A for hxed tc also tends to reduce the typi- 
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Figure 7. (a) The mean Lyapunov exponent (x) (upper two 
curves) and the associated dispersion (lower two curves) com¬ 
puted for ensembles of 1600 orbits evolved in an undamped Plum¬ 
mer potential with amplitude m = 0.5, ojq = 3.5, and variable 
autocorrelation time tc- Diamonds correspond to = 0.35, tri¬ 
angles to = 2.8. (b) The same for mo = 0.1. (c) mo = 0.05. 


cal size of a finite time X- This, however, is not generic. If, 
instead, the same computations are repeated for mo = 0.05, 
one finds that the choice of A which leads to the largest x’s 
actually varies with tc. 

The effects of varying A and tc on the mean (x) and the 
dispersion associated with these ensembles is exhibited 
in Figure 7. 


4.2 Damped oscillations 

Allowing for damped oscillations of the form given by eq. (5) 
or (6) has a comparatively minimal effect. Not surprisingly, 
the size of the mean Lyapunov exponent (x) decreases as 
time elapses and there can be a gradual decrease in the rel¬ 
ative fraction / of orbits exhibiting exponential sensitivity. 
However, the basic phenomenon of transient chaos persists 




Figure 8. (a) The four lower curves exhibit the mean value (x) as 
a function of ui for chaotic orbits subjected to damped oscillations 
with mo = 0.5, p = 2, and to for the intervals (from top to 
bottom) 100 < t < 200, 200 < t < 300, 300 < t < 400, and 
400 < t < 500. The top curve exhibits (x) for undamped (p = 0) 
oscillations for the intervals 0 < t < 512. (b) The same data 
plotted as a function of logjQ uj. 


until the perturbation has damped almost completely away. 

This is, e.g., evident from Figure 8, which was again 
generated from the same 1600 initial conditions used to gen¬ 
erate Figures 1-5, now allowing for a perturbation of the 
form given by eq. (6) with mo = 0.5, to = 100, and p = 2, 
again allowing for variable lo. Here the four lower curves 
(from top to bottom) represent mean values of the finite 
time Lyapunov exponent (x) generated separately for the 
intervals 100 < t < 200, 200 < t < 300, 300 < t < 400, 
and 400 < t < 500. The top curve, reproduced from Fig. 2, 
represents (x) for 0 < t < 512 for the same orbits evolved 
in the presence of undamped {p = 0) oscillations with the 
same initial mo. It is evident that, as one would have ex¬ 
pected, the degree of exponential sensitivity decreases with 
time and that, by the end of the integration, the orbits are 
behaving in a nearly regular fashion. The fact that the peak 
frequency lo drifts towards lower values at later times is con¬ 
sistent with the fact, evident from Fig. 2, that for undamped 
oscillations, lower amplitudes mo correlate with lower values 
of UJ. 

That it may be possible to achieve efficient chaotic 
phase mixing in an oscillating galactic potential while still 
relaxing towards a nearly integrable state within 10 — 20tD 
or so is illustrated visually in Figures 10 - 12, which exhibit 
the behaviour of the same localised ensemble of 1600 initial 
conditions, again subjected to damped oscillations of the 
form given by eq. (6) with mo = 0.5, p = 2, and to — 100 
for three different choices of frequency uj. Figure 9 exhibits 
the results of a corresponding evolution in the presence of 
a nonoscillatory perturbation of the form given by eq. (8). 


© 2000 RAS, MNRAS 000, 1-?? 












Transient chaos and resonant phase mixing in violent relaxation 9 


4 

2 

> 0 

-2 

-4 

( o ) 

4 

2 

^ 0 

-2 

-4 

( b ) 


CN 

O 

CN 

1 


4-20 2 ^ 


X 


X 

4 


4 


2 


2 






> 0 


^ 0 


-2 


-2 



( c ) 


( d ) 

-4 


-4 



CN 

O 

CN 

1 


4-20 2 ^ 


X 


X 

4 


4 


2 


2 


> 0 

r 

^ 0 


-2 


-2 



( e ) 


(0 

-4 


-4 



-4 -2 0 2 4 -4 -2 0 2 4 

X X 


Figure 9. The x and y coordinates of an initially localised en¬ 
semble of orbits evolved in a Plummer potential subjected to a 
nonoscillatory perturbation of the form given by eq. (7) with 
mo = 0.5 and p = 2. (a) t = 0. (b) t = 16. (c) t = 32. (d) 
t = 64. (e) t = 128. (f) t = 256. 


In each case, the six panels exhibit the x and y coordinates 
of each of the orbits at times t = 0, 16, 32, 64, 128, and 
256, the last corresponding to an interval of approximately 
12.8tD. 

It is evident visually from Figure 9 that, in the ab¬ 
sence of oscillations, mixing is comparatively inefficient. In¬ 
deed, the evolution in that Figure is qualitatively identical 
to examples of regular phase mixing in time-independent 
potentials (c/. Figure 2 in Kandrup 1999). As asserted in 
Section 3, a non-oscillatory perturbation of the integrable 
Plummer potential leads to little if any transient chaos and, 
as such, no evidence for chaotic phase mixing. The remain¬ 
ing three Figures, which incorporate a systematic pulsation, 
all yield phase mixing that is substantially more efficient. 
Figure 10, which represents orbits that have been pulsed 
with LO = 0.035, a frequency near the lower edge of the res¬ 
onance, clearly exhibits more robust mixing, although the 
mixing is still considerably less efficient than what is ob¬ 
served for strongly chaotic flows in time-independent poten¬ 
tials (c/. Figure 1 in Kandrup 1999). Indeed, it is difficult to 
determine unambiguously whether the orbits used to gener¬ 
ate this Figure genuinely exhibit significant transient chaos. 

By contrast. Figure 11, which was generated for w = 
0.70, exhibits precisely the sort of behaviour which one 
has come to associate with chaotic phase mixing in time- 
independent potentials. For the first two dynamical times, 
the localised ensemble of orbits, which started with a phase 
space size < 10“^, still remains confined. By t = 64, how- 
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Figure 10. The same as the preceding Figure, now allowing for 
an oscillatory perturbation of the form given by eq. (6) with Mq = 
0.5, p = 2, and u! = 0.035. 

ever, corresponding to an interval « 3.2t_D, the orbits have 
begun to spread significantly; and, hy t = 128, correspond¬ 
ing to roughly 6.4t_D, the ensemble has dispersed to fill a 
large portion of the available phase space. Nevertheless, de¬ 
spite this efficient mixing, the potential is damping rapidly 
towards a near integrable state. By a time t = 256, the 
overall amplitude of the perturbation mo/(l -I- t/to)^, has 
decreased from mo = 0.5 to mo ~ 0.04, this corresponding 
to a much more nearly time-independent system. 

Figure 12 exhibits an example of ‘incomplete’ chaotic 
phase mixing for the case of a frequency w = 3.5. Here again 
one sees clear evidence of the dispersive behaviour indicative 
of chaotic phase mixing; but, since one is closer to the up¬ 
per edge of the resonance, the efficient mixing has de facto 
turned off before the ensemble could fill all the accessible 
phase space. The dearth of orbits near a: = y = 0 is exactly 
what one would expect in a time-independent spherical po¬ 
tential for orbits with an appreciable angular momentum. 


5 DISCUSSION 

The numerical computations described herein lead to several 
seemingly unambiguous conclusions: 

1. Subjecting orbits to a possibly damped oscillatory time- 
dependence can give rise to substantial amounts of transient 
chaos, even if the initial and final form of the potential is 
completely integrable. 

2. This transient chaos appears to arise from a resonant cou¬ 
pling which, at least for large amplitudes, is very broad and. 
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Figure 11. The same as the preceding Figure, for lo = 0.70. 



Figure 12. The same as the preceding Figure, for lj = 3.50. 


as such, might be expected to be present in many physical 
systems. 

3. This transient chaos can drive chaotic phase mixing which, 
in the context of a fully self-consistent evolution, might be 
expected to play an important role in violent relaxation, 
e.g., during the formation of galaxies and galaxy halos and 
in mergers/close encounters of galaxies. 

Three extensions of this work are underway: The first 
involves a consideration of additional examples with the aim 
of determining the extent to which this transient chaos is 
generic. The examples described here suggest that transient 
chaos can arise whenever the perturbation has substantial 
power at frequencies that are comparable to within an order 
of magnitude or so to the natural frequencies associated with 
the orbits. However, this needs to be checked in the context 
of other, more realistic models. 

A second extension involves efforts to better understand 
the origins of this chaos using geometric arguments due orig¬ 
inally to Pettini (1993). Recently, Kandrup, Sideris, & Bohn 
(2002) showed that thermodynamic arguments proposed 
originally to predict the size of the largest Lyapunov expo¬ 
nent in high-dimensional time-independent Hamiltonian sys¬ 
tems (Casetti, dementi, & Pettini 1996) also work surpris¬ 
ingly well in lower-dimensional {e.g., two- and three-degree- 
of-freedom) systems. The aim is to adapt that work to the 
case of Hamiltonian systems admitting a time-dependence of 
the form that one might expect to encounter in the context 
of violent relaxation. 

The third extension involves searching for evidence of 
transient chaos in the context of fully self-consistent numer¬ 
ical simulations, both for self-gravitating systems and for 
non-neutral plasmas, e.g., charged particle beams. Recent 
simulations {e.g., Kishek, Bohn, Haber, O’Shea, Reiser, & 
Kandrup 2001) have shown that fully self-consistent simu¬ 
lations of beams can exhibit evidence of chaotic phase mix¬ 
ing. The obvious issue is whether, as one might expect, the 
degree of chaotic mixing observed in such a beam or in a 
simulation of violent relaxation in a self-gravitating system 
correlates with the degree to which the bulk potential admits 
a significant time-dependent oscillatory component. 

Plummer spheres do not constitute especially realistic 
models of early-type galaxies: aside from the idealisation of 
spherical symmetry, real galaxies often have a central cusp 
and their densities fall off much more slowly at large radii. 
Despite this, however, it is interesting to estimate more care¬ 
fully the ‘real’ values that H and lo might assume if nature 
did somehow construct Plummer spheres. 

To estimate a reasonable value for the pulsation fre¬ 
quency one can proceed as follows: In the dimensionless 
units used in this paper, the half-mass radius rn inside 
of which half the mass of the galaxy is contained satisfies 
Th, ~ 1.3. However, r{E), the mean radius associated with 
a uniform sampling of an S = constant hypersurface, sat¬ 
isfies r{E) = Th for E « 0.4. One might thus argue that a 
‘typical’ frequency associated with the system as a whole 
would correspond to a characteristic frequency associated of 
orbits with E ~ 0.4 which, by analogy with Figure 5, can be 
computed as Hm ~ 0.4. Alternatively, given that the den¬ 
sity p = (3/47r)(l -I- it is easily seen that the half 

mass density ph = p{rh) defines a characteristic frequency 
^Ih = V^'^Gph ~ 0.50. One might, therefore, anticipate that 
realistic bulk oscillations associated with the system would 
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be characterised by frequencies 0.4 ^ u} ^ 1.0, the larger val¬ 
ues corresponding to ‘higher order modes.’ Alternatively, an 
examination of Fourier spectra for representative orbit en¬ 
sembles at different energies shows that the peak frequency 
ranges from Qrn 1.0 for E —1.0 to, e.g., Q.m ~ 0.05 for 
E = -0.1. 

To the extent that the resonance described in this paper 
is representative of the Plummer sphere, it would thus seem 
likely that all but the very highest energy orbits would in 
fact be able to resonate with a large scale ‘bulk’ oscillation. 
One might also be tempted to conclude that the possibility 
of resonances for uj ^ is largely unimportant physically. 
This, however, is not necessarily the case. One anticipates 
generically that, as the system ‘relaxes’, power will cascade 
from larger scales and lower frequencies to shorter scales and 
higher frequencies. To the extent, however, that the natu¬ 
ral frequency of the perturbation is more important than its 
detailed form, one might expect that such higher frequency 
oscillations, once triggered, could play an important role in 
continuing the process of violent relaxation on shorter scales. 
This is, e.g., consistent with the fact (Kandrup et al 2003) 
that shorter scale perturbations associated with supermas- 
sive black hole binaries orbiting near the center of a galaxy 
can trigger qualitatively similar effects in terms of resonant 
phase mixing and energy and mass transport. 
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